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The nature of light-induced structures in attosecond transient absorption spectroscopy of molec¬ 
ular systems is investigated theoretically. It is shown how nuclear dynamics affect these structures. 
We find that a theoretical three-surface model captures the main characteristics in the calculated 
spectra. Based on this model, nuclear dynamics is divided into different categories, each category 
having unique signatures in the absorption spectra. Finally, we discuss the possibility for experi¬ 
mental observation of light-induced structures in molecules. 
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I. INTRODUCTION 


Since 2010 attosecond transient absorption spec¬ 
troscopy (ATAS) has been used to study electronic dy¬ 
namics in atoms on the attosecond time scale urn- Dur¬ 
ing propagation through a target region affected by an 
infrared (IR) or a near-infrared (NIR) pulse, the spectral 
profile of the attosecond pulse is modified. After propa¬ 
gation the intensity spectrum is recorded to measure the 
accumulated modification. This spectrum depends on 
the relative phase between the attosecond pulse and the 
NIR/IR pulse and therefore on the time delay between 
the two pulses. Indirect temporal information about the 
response of the target to the NIR/IR field is therefore 
obtained from the delay-dependent intensity spectrum. 
The spectral profile of the attosecond pulse is in the ex¬ 
treme ultraviolet (XUV) regime. The intensity of the 
XUV pulse is low, and only processes involving a sin¬ 
gle XUV photon are non-negligible. The intensity of the 
NIR/IR pulse is in a regime where few-photon processes 
occur, but the intensity is too low to give rise to a signal 
in the XUV regime of interest. 

Of particular interest in ATAS is the time interval 
where the XUV and the NIR/IR pulses overlap. In this 
region a number of spectral features arise in the spec¬ 
trum as a consequence of non-linear coupling between the 
system and the two fields [ME]. Here we focus on light- 
induced structures (LISs). LISs are signatures of non¬ 
linear processes transferring population from the ground 
state to a state of the same parity (a dark state) by 
absorbing an XUV photon and absorbing or emitting 
a NIR/IR photon. In atomic systems LISs are there¬ 
fore seen in ATAS spectra as pairs with each individual 
absorption feature located one NIR/IR photon energy 
above (emission of NIR photon) or below (absorption of 
NIR photon) the dark state energies. The spectral width 
of an atomic LIS is inversely proportional to the duration 
of the NIR/IR pulse and since this pulse is often almost 
as long as the natural damping of the system, LISs ap¬ 
pear as relatively narrow features in the ATAS spectrum 
O [9] . In Ref. [TO] it was shown that molecular ATAS 
spectra are to a large extend determined by molecular 


vibrations, and the results indicated that LISs at well- 
defined energies are not a general molecular property. It 
is therefore of interest to make a careful investigation of 
how the LISs are affected by nuclear motion. It turns 
out that molecular LISs are extremely sensitive to the 
the type of nuclear dynamics supported by the electronic 
dark state Born-Oppenheimer (BO) curve. LISs there¬ 
fore lend themselves to a study of nuclear dynamics in a 
new setting and they form a possibility to study the com¬ 
bined nuclear and electronic dynamics in the otherwise 
hidden dark states. The main purpose of this work is 
to give a more complete description of LISs in molecules, 
and show how LISs can be used to classify molecular pro¬ 
cesses. 

The paper is organized as follows. Section [H] intro¬ 
duces the theoretical model used to obtain the ATAS 
spectra, the systems, and parameters used in the calcu¬ 
lations. Further, we introduce a three-curve model to iso¬ 
late the main physical mechanics of LISs in molecules. In 
Sec. Ill we first show LISs in the ATAS spectrum of H 2 ’^ 
for fixed and moving nuclei. In H 2 ^ all excited BO curves 
are dissociative. To investigate the effect of bound nu¬ 
clear motion on ATAS spectra in general, and on LISs in 
particular, we consider classes of generic potential energy 
curves. The resulting ATAS spectra are simple enough 
that the characteristics of LISs can be identified clearly. 
Section |IV| concludes the paper and relates the results to 
a possible experimental observation of LISs in molecules. 
Atomic units are used throughout. 


II. MODELS 

To calculate the delay-dependent ATAS spectrum we 
use the method of Ref. m (see also Ref. [10] where this 
method was related to other approaches). Assuming that 
effects accumulated in propagation though the target re¬ 
gion can be captured in a single-system response model 
[12], the function 

= 2Re[£’*„(w,T)£’gen(w,T)] , (1) 

gives the modification of the intensity spectrum at angu¬ 
lar frequency uu and time delay r. In Eq. 0- Ein(w,r) is 
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the Fourier transform of the unmodified, incoming field 
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where n is the target density, c is the speed of light, 
and ((i(t, r)) is the expectation value of the single-system 
dipole moment, found from quantum mechanical calcu¬ 
lations. 

We consider fields linearly polarized along the 2 ;- 
direction, and propagating in the positive x-direction. 
The 2 ;-component of the electric field E{t) = -dtA{t) 
used to model the XUV and NIR/IR pulses is defined 
through the vector potential 


A{t) = Aq exp 


{t-tcf 


COS [uj{t - tc)] , 


T = NcTc = Nc — , 

UJ 


(4) 


where Aq = and Eq relates to the intensity as 

I = with one atomic unit of intensity equal to 

3.51 X 10^^ W/cm^. In Eq. Q tc is the center of the 
pulse. The full width at half maximum ^WHM) dura¬ 
tion TpwHM is related to T used in Eq. 0 by TpwHM = 
To model H 2 '' we consider a simplified one¬ 
dimensional (ID), soft Coulomb potential model where 
the molecular axis is aligned with the linearly polarized 
laser pulse. Within this model, electronic and nuclear 
degrees of freedom are treated exactly. Eor more infor¬ 
mation see Refs. [niiii]. To mimic experimental results 
we damp the dipole moment using a window function of 
a duration of 42 fs m- 

In Ref. [To] it was shown that the dynamics in H 2 ’^ of 
importance for ATAS are well-described by an expansion 
in a few electronic states. In this work we therefore ex¬ 
pand the wave functions in the lowest N electronic BO 
states, (j)j{z]R)\ 


N 

j = l 

where Gj{RA) is the nuclear wave packet corresponding 
to the electronic state (pj{z; R), R is the internuclear sep¬ 
aration, and z denotes the electronic coordinate. Eor the 
analysis of nuclear dynamics it is sometimes useful to fur¬ 
ther expand the nuclear wave packets Gi{RA) in vibra¬ 
tional eigenstates Xi,k{E) (discrete states) and XEi{R) 


(nuclear continuum states) corresponding to the i’th elec¬ 
tronic state 

k 

+ f dEi CEMxEAR)e-'^^' 

= ^dEi^k Ci,k{t)xi,k(R)e~''^'''"\ ( 6 ) 

where Ei^^ and Ei are the discrete and the continuous 
eigenenergies of the vibrational states of the Fth elec¬ 
tronic state. Note that all vibrational states correspond¬ 
ing to a specific electronic state can be nuclear-continuum 
states (see Ref. uni for details). 

LISs found in ATAS spectra are a consequence of the 
NIR/IR coupling between an excited state of parity op¬ 
posite from the ground state (bright state) and an ex¬ 
cited state with the same parity as the ground state; the 
dark state. The latter state is refereed to as a dark state 
since the transition from the ground state to this state 
is dipole forbidden and there are therefore no spectral 
features in the dipole response corresponding to the en¬ 
ergy difference between dark states and the ground state. 
Consequently these states are not observed in one-photon 
spectroscopy. In Ref. [5] it was shown that the main char¬ 
acteristics of LISs in atoms are often captured in calcu¬ 
lations using a three-state model. Even though LISs and 
the effect of nuclear motion on LISs was not the main fo¬ 
cus of Ref. nni, results of that paper indicated that LISs 
in molecules are also often well-described in a model only 
including three electronic states (three-curve model). In 
the three-curve model the ground state (j)g{z^R)^ the 
bright state (j)c{z^R)^ and the dark state (j)d{z^R) [see 
Eq. <©] are included. Within this model a general state 
can therefore be written as 


-^{z^RR) =Gg{R,t)(f>g{z-, R) + Gd{R,t)(pd{z\R) 

+ Ge{R,t')(pe{z‘, R) 

=Xg,o{R)4>g{z', R) 

+ ^dEd,iCd,iit)xdARy~'^"-‘^'~"Ud{z; R) 

+ dEe,kCe,k{t)Xe,k{Ry^^''‘^'~^Ue{z-, R). 

(7) 

In obtaining the second part of Eq. 0 it is used that 
the ground state to a good approximation remains un¬ 
changed during propagation. To analyze LISs in molec¬ 
ular ATAS spectra we also use the perturbative limit of 
the three-curve model: Assuming that the electric field 
is monochromatic [E{t) = Eq sin(u;t)] we obtain that the 
first-order change in the expansion coefficient Ce^k{t)^ for 
the Eth vibrational state of the bright electronic state, 
is given by 



EodtARo) 

2i 



” ^i(Ee,k-Ed,i+uj)x(t-T) 

Ee,k - EdA + ^ 


I 











3 




Ee,k - Ed^i - LU 


{Xe,k\Xd,l)cd,l{T )■ 

( 8 ) 


The expectation value of the dipole moment in the three- 
curve model of Eq. 0 is well approximated by 


(d(t, t)) ~2Re 



{Xg,o\Xe,k) Ce,/c(^) 
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( 9 ) 



To obtain Eq. 0 it is assumed that the electronic dipole 
moment function 

dfy{R) = J dz4>*{z-R)d<pj{z-,R) ( 10 ) 

is a slowly varying function of R and therefore can be 
evaluated at Rq = f R\Gg{RR = -oo)\‘^dR. This is usu¬ 
ally a good approximation for ATAS calculations (see 
Ref. [iQ]). Inserting Eq. ^ into Eq. ^ shows that the 
energies disappear in the final expression for the 

time-dependent dipole moment, leaving terms oscillating 
with frequencies ± cjc- The distribution and popu¬ 
lation of vibrational states corresponding to the bright 
electronic state (j)e{z^R) does therefore not affect the de¬ 
phasing of the dipole moment responsible for LISs. The 
distribution and population of vibrational states corre¬ 
sponding to the dark electronic state (j)d{z^R) are, on 
the other hand, crucial for the nature of LISs. The more 
dark vibrational states that are populated, the faster the 
resulting dipole moment will dephase. Even in the situa¬ 
tion where a finite number of discrete vibrational states 
are populated, and one would expect a partial revival 
of the time-dependent dipole moment, the revival is of¬ 
ten absent due to the finite duration of the NIR pulse. 
We therefore expect LISs with well-defined energies to be 
suppressed in the molecular ATAS spectrum when more 
than a few vibrational dark states are populated. The 
results in the next section verify this explanation. 


III. RESULTS 


In this section we present ATAS spectra calculated us¬ 
ing Eq. 0- To find the time-dependent dipole moment 
we propagate the nuclear wave packets of Eq. 0 in ac¬ 
cordance with the time-dependent Schrodinger equation 
(see Ref. m for details). Eirst we calculate the ATAS 
spectra for our H 2 ’^ model with fixed and moving nu¬ 


clei, then in Sec. IIIB we investigate how LISs depend 
on the relative arrangement and shapes of the electronic 
BO curves. 


A. H2^ model 

In Eig. we show the lowest six BO potential curves 
of our ID model of H 2 ^. In the model we use a soft- 


R (a.n.) 

Figure 1. BO curves of our H 2 ^ model. The thick, full 
line shows the ground state, while the narrower full lines are 
bright excited states and the dashed lines are excited dark 
states (see text). The vertical gray line indicates the value of 
Ro = f R\Gg(RR = -oo)\^dR = 2m. 


ened Coulomb potential designed to reproduce the exact 
three-dimensional IscFg BO curve [131 US]- The figure 
shows that a nuclear wave packet located on any excited 
BO curve dissociates. In Eig. the ATAS spectra are 
shown for the atomic-like situation with the nuclei fixed 
at R = Rq (left panel) and for the molecular case includ¬ 
ing nuclear motion (right panel). We immediately see 
from Eig.j^that including nuclear motion in the descrip¬ 
tion completely changes the ATAS spectra from atomic- 
like spectra characterized by narrow absorption lines to 
spectra dominated by very broad features. In the fol¬ 
lowing we will see that such broadening is characteristic, 
not only for spectral features due to dipole allowed tran¬ 
sitions, but also for LISs. We investigate how selected 
electronic states affect the spectrum by eliminating them 
from the basis set of the calculation. In this way the 
physical processes responsible for different spectral fea¬ 
tures are identified. 

The ATAS spectra of Eig.j^a) are calculated using the 
expansion 0 including the ground state and the five low¬ 
est excited electronic states. The energies of the four elec¬ 
tronic states in the energy range of relevance are shown 
in the middle column of Eig. [^a), where energies corre¬ 
sponding to bright and dark states are indicated by full or 
dashed lines, respectively. In Eigs.[^b)-(e) some of these 
lines are shown in light-gray coloring, indicating that the 
electronic state of that energy has been omitted in the 
calculations of the spectra. All energies are given with 
respect to the ground state energy. The fixed nuclei spec¬ 
trum [Eig.|^a), left] shows many similarities with atomic 
spectra. Eor large positive delays of the T = 330 as XUV 
pulse [see Eq. Q] with respect to the NIR pulse the spec¬ 
trum is characterized by narrow absorption features at 
the field-free energies of bright excited states at 21.48 eV 
and 25.29 eV (full lines in the middle column). The LISs 
in the fixed nuclei spectrum, located one NIR-photon en¬ 
ergy above/below the dark state of energy 17.65 eV, are 
encircled by black cigar shapes and are characterized by 
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droplet-like shapes, a periodic modulation with a period 
of half the period of the NIR cycle, and a width inversely 
proportional to the duration of the NIR pulse. These are 
all characteristics known from LISs in atomic systems 
We use a Tnir=14 fs NIR pulse [see Eq. @] and 
as a consequence our LISs have widths in the sub-eV do¬ 
main. In the moving-nuclei spectrum the situation is very 
different [Fig.j^a), right]. All narrow absorption features 
are gone, and the two absorption lines corresponding to 
transitions to bright states (at 21.48 eV and 25.29 eV in 
the fixed-nuclei spectrum) have merged into a single very 
broad feature. More important for the current discussion 
is that the LISs at the relatively well-defined energies of 
~ 16 and ~ 19 eV appear to be absent in the moving-nuclei 
spectrum. The only visible NIR-field modulation of the 
moving-nuclei spectrum in the right panel of Fig. I^a) is 
a periodic modulation with a period of half a NIR-field 
cycle from -5 fs to 5 fs, extending over the entire energy 
range of the figure. Since this modulation seems to be 
centered around the broad absorption line at ~ 21-24 
eV it is natural to assume that this modulation is a con¬ 
sequence of ’which-way interference’ [8] between the two 
bright states [at 21.48 eV and 25.29 eV, indicated by full 
lines in the middle column in Fig. m- In Fig.^b) we 
calculate the ATAS spectrum again, but this time with¬ 
out the highest excited bright state (now a gray line at 
25.29 eV in the middle column). Since only one bright 
state is left in the frequency domain of interest, we ex¬ 
pect all modulations caused by which-way interference 
to be eliminated. From the moving-nuclei spectrum of 
Fig-Hb), however, we see that the spectrum is remark¬ 
ably unaffected by the removal of the bright state. For 
the fixed-nuclei spectrum of Fig. [^b) we see that the 
LISs at ~ 16 eV and ~ 19 eV are also unaffected. Since 
the interference between bright states apparently can not 
explain the modulation of the moving-nuclei ATAS spec¬ 
trum we are led to expect that this modulation is caused 
by LISs. To validate this expectation we calculate the 
ATAS spectrum without the highest excited dark state 
[at 23.87 eV, gray dashed line of Fig. |^c)] and without 
the second highest excited dark state [at 17.65 eV, gray 
dashed line of Fig. [^d)]. Without the dark states there 
will not be any LISs, and from the fixed-nuclei spectrum 
of Fig. [^d) we see that the LISs at ~ 16 eV and ~ 19 
eV are indeed missing when the 17.65 eV dark state is 
removed from the basis set of the calculation. In the 
moving-nuclei spectra we see that removal of the highest 
excited dark state reduces the modulation in the high 
energy part of the spectrum [Fig. [^c)], and removal of 
the lower excited dark state reduces the modulation in 
the low energy part of the spectrum [Fig. [^d)]. Finally 
in Fig. I^e) we remove both dark states from the cal¬ 
culation. In the corresponding fixed and moving-nuclei 
spectra the effect of the NIR pulse is almost gone. It 
is not surprising that Fig. |^e) is almost independent of 
the NIR pulse since all dark states have been removed. 
When Fig. I^e) is hold together with Figs, [^b)-(d) it is, 
however, clear that the NIR-modulations in the moving 
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Figure 2. (Color online) ATAS spectra S{uj,r) [Eq. Q] 
of H 2 ^ with fixed (left column) and moving (right column) 
nuclei, (a) includes the six lowest electronic states in the 
calculation. The energy evaluated at Rq of the four states in 
the energy range of interest are shown in the middle column. 
In the middle column full lines represent bright states and 
dashed lines represent dark states (see text), (b)-(e) excludes 
the electronic states indicated by light-gray coloring in the 
middle column. Note that there are separate color bars for 
the left and right columns. The color bars are designed to 
highlight the LISs which are also encircled by black cigars 
in the fixed nuclei case. Pulse parameters: Txuv = 330 as, 
Tnir = 14.0 fs, Axuv = 50 nm, Anir = 700 nm, /xuv = 5 x 10^ 
W/cm^ /nir = 2 X 10^^ W/cm^ 
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nuclei spectrum of Fig.|^a) are primarily resulting from 
LISs. 

Returning to Fig. |^c) we observe that the removal of 
the highest excited dark state causes a drastic reduction 
in the Tnir/ 2 oscillations in the absorption lines at 21.48 
eV and 25.29 eV. These oscillations are a result of which- 
way interference, and we therefore conclude that the dark 
state is essential in the coupling of the two bright states 
which is also anticipated from, e.g., second-order pertur¬ 
bation theory. Additionally it is seen from Fig. [^c) that 
the absorption feature at ~ 23 eV stretching from —10 
to ~ 5 fs [in Figs.j^a), (b) and (d)] disappears from the 
spectrum. This absorption feature is therefore also a LIS. 

We have argued that the periodic modulation of the 
moving-nuclei spectrum of H 2 ’^ [Fig.j^a), right], stretch¬ 
ing over a broad energy range, are in fact LISs. To ex¬ 
plain why these LISs do not look anything like what we 
see in atomic systems and as in H 2 ’^ with fixed R we re¬ 
turn to Eqs. and From these equations we were 
able to conclude that LISs would be broadened if more 
dark vibrational states were populated since this leads to 
a dephasing of the part of the time-dependent dipole mo¬ 
ment responsible for LISs. Based on the Frank-Condon 
principle we expect that the nuclear wave function, ex¬ 
cited to a dark state BO curve, will initially be local¬ 
ized around R = Rq, but since all dark states in H 2 ’^ are 
dissociative, all the corresponding vibrational states are 
nuclear-continuum states. An expansion of the nuclear 
wave packet in vibrational dark states will therefore in 
principle require an infinite number of dark vibrational 
states. As a result we expect that the corresponding 
dipole signal will be damped on a very short time scale 
leading to extremely broad spectral features. How broad 
the LISs are depends on the distribution of vibrational 
dark states populated by the combined XUV and NIR 
pulses and therefore on the slope of the corresponding 
BO curve. 


B. Illustrative examples using the three-curve 
model 

In many situations LISs can be accurately described 
using a three-curve model containing the electronic 
ground state, one electronic dark state and a single elec¬ 
tronic bright state [see Eq. 0]. We take the ground state 
BO curve to be binding. Eor the excited curves, we con¬ 
sider different shapes. Our goal is to investigate how the 
shapes of the two excited BO curves in the three-curve 
model affect the LISs. We divide the excited BO curves 
into two categories: binding and dissociative. This di¬ 
vision gives us four different characteristic arrangements 
of BO curves as shown schematically in Eig. The up¬ 
per left set of BO curves in Eig. represents the case 
where both excited BO curves are dissociative. This case 
is similar to H 2 ’^ discussed in the previous section. We 
have confirmed that the three-curve model is capable of 
reproducing the LISs in the ATAS spectra of Eig. cor- 



Figure 3. (Color online) Illustration of possible arrangements 
of excited BO curves in molecules. Full lines represent bright 
states with respect to coupling from the ground state curve, 
while dashed lines represent dark states. 


responding to the included dark state, for both fixed and 
moving-nuclei calculations. To obtain the other arrange¬ 
ments of BO curves shown in Fig. we replace one or 
two of the excited BO curves Cd{R) (dashed, dark state 
curve) and Ce{R) (full, bright state curve) of H 2 ^ with 
the curves 


= Cg{R - Ra/e) - e,(i?o) + (11) 


where Cg{R) is the ground state curve. For R^je = 0 
the curve Cd/e{R) is parallel to the ground state and the 
energy difference to the ground state curve is as for the 
original curve Cdfe{R) evaluated at Rq. In the numerical 
implementation of the model we have used the dipole ma¬ 
trix elements from the H 2 ^ case. We have found that the 
main characteristics of LISs in the spectra of the three- 
curve model are not largely affected by the specific choice 
of dipole matrix element values, when reasonable values 
are chosen. 

The upper right arrangement of BO curves in Fig. 
represent the case where both excited curves are bind¬ 
ing. In Fig. we show the corresponding ATAS spectra 
for different values of R^ and Re? be., for different rel¬ 
ative positions of the bright and dark state BO energy 
curves. In Fig. ga) {Rd = Re = 0) we see a narrow ab¬ 
sorption line at ~ 22 eV for the dipole allowed transition 
and a narrow single-line LIS at ~ 19 eV. This atomic-like 
behavior is expected since all three BO curves are par¬ 
allel. The vibrational ground state is therefore the same 
in all three curves, and from the Franck-Condon princi¬ 
ple nuclear dynamics does not play any important role, 
since excited vibrational states are never populated. In 
Fig.|^b) the minimum of the bright state BO curve is at 
larger R [Eq. (II) with Re = 0.3]. When the lowest vibra¬ 
tional state of the electronic ground state is expanded in 
the (discrete) vibrational states of the bright state curve 
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Figure 4. (Color online) ATAS spectra of a model with two 
binding excited BO curves [see upper right part of Fig. [^. 
The spectra are calculated for the different values of Rd and 
Re [see Eq. (11)] Rd = Re = 0.0 (a), Rd = 0.0 and Re = 0.3 
(b), Rd = 0.3 and Re = 0.0 (c), Rd = 0.3 and Re = 0.3 (d), 
Rd = 0.5 and Re = 0.0 (e), and Rd = 0.5 and Re = 0.3 (f). 
Special attention should be given to the spectral feature in (a) 
and (b) at ~ 19 eV. The insets show illustrations of the two 
excited BO curves included in the calculations. In the insets, 
the vertical gray line indicate R = Rq. Pulse parameters are 
as in Fig. 


several of the expansion coefficients will be non-zero. In 
the spectrum we therefore see several absorption lines 
corresponding to the energy of the individual vibrational 
states. The LIS at ~ 19 eV in Fig. [^b), however, are 
remarkably unaffected by the translation of the bright 
state curve. This is in agreement with the conclusion 
in Sec. [Ilj that the LISs are unaffected by the popu¬ 
lation and distribution of vibrational bright states. In 
Figs. [^c) and m the minimum of the dark state BO 
curve is moved to larger R in the calculations {Rd = 0.3). 



Figure 5. (Color online) ATAS spectra of a model with a sin¬ 
gle binding excited bright BO curve and a single dissociative 
dark BO curve [see lower left part of Fig. [3. In the calcula¬ 
tions Re of Eq. (11) is chosen as = 0 (a) or Re = 0.3 (b). 
The insets show illustrations of the two excited BO curves 
included in the calculations. In the insets, the vertical gray 
lines indicate R = Rq. Pulse parameters are as in Eig. 


As a result several vibrational dark states are populated 
and we observe a broadening of LISs in the spectrum. A 
careful examination of Figs. |^c) andj^d) shows that the 
LIS at ~ 19 eV has a horizontal few-line structure (not 
clearly visible at the scale of the figure) corresponding to 
the energies of the different populated vibrational dark 
states. We therefore refer to this type of spectral fea¬ 
tures as few-line molecular LISs. In Figs. I^e) and if) 
the spectra for calculations with Rd = 0.5 are shown, and 
we see even broader absorption lines without horizontal 
line structure (defuse LISs). In the latter case the slope 
of the dark BO curve evaluated at Rq is steeper than for 
the Rd = 0.3 case [Figs. ic) and[^d)] and it is therefore 
expected that more vibrational states become populated 
by the combined action of the XUV and NIR laser pulses. 
The absorption features in Figs, [^a)-(f) corresponding 
to dipole allowed transitions seems to be unaffected by 
the dark state. However, we know from the discussion 
of Fig. [^that dark states play an important role in the 
transitions between bright states. In general, the shapes 
of dark state BO curves therefore affect the absorption 
lines corresponding to dipole allowed transitions. 

We now investigate the case of a bound bright BO 
curve and a dissociative dark BO curve schematically 
drawn in the lower left part of Fig. The correspond¬ 
ing ATAS spectra in Fig. show that absorption lines 
resulting from dipole allowed transitions are again very 
similar to what we saw in Fig. The LISs, however, are 
very defuse, stretching over a broad energy region. This 
is similar to what we saw in the movingmuclei spectra 
of H 2 ’^ shown in the right panel of Fig. ^a). The fact 
that the LISs of Fig. are similar to LISs in the moving 
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Figure 6. (Color online) ATAS spectra of a model with a sin¬ 
gle dissociative excited bright BO curve and a single binding 
dark BO curve [see lower right part of Fig. [3. In the calcu¬ 
lations Rd of Eq. 0 is chosen as = 0 (ajor = 0.3 (b). 
The insets show illustrations of the two excited BO curves 
included in the calculations. In the insets, the vertical gray 
lines indicate R = Rq. Pulse parameters are as in Fig. 


nuclei spectrum of H 2 ’^ is another example that often the 
characteristics of LISs do not depend on the shape of the 
bright state curve. 

In Fig. we show the spectra of a molecule with a 
bound dark state curve and a dissociative bright state 
curve as sketched in the lower right corner of Fig. In 
Fig. [^a) we see that there are several absorption fea¬ 
tures above the main LIS at ~ 19 eV [also seen in, e.g., 
Fig. I^a)]. Since the ground state and dark state curves 
are parallel, this structure indicates that the LISs are in 
fact affected by the bright state curve in this case. For 
the calculations carried out to produce Fig.j^a) the slope 
of the dark state BO curve evaluated at Rq is so steep 
that the nuclear wave function Ge(R,t) slide down the 
curve on a timescale comparable to the duration of the 
XUV pulse. As a result several vibrational dark states 
are populated by the combined action of the XUV and 
NIR fields even though the ground state and dark state 
curves are parallel. In spectra for models where a less 
steep dissociative bright state is used, the LISs are again 
unaffected by the shape of the bright state curve. For 
the calculations producing Fig.j^b) the minimum of the 
dark state curve has been moved to larger R-values, but 
the spectra of Figs, [^a) and |^b) are very similar. In 
contrast to earlier cases it is difficult to say much about 
the position and shape of the dark BO curve from LISs 
when the bright state curve is very steep compared to 
the dark state curve. The effect described above can be 
minimized by using a shorter XUV pulse. 


IV. CONCLUSION AND DISCUSSION 


We have investigated LISs in four different arrange¬ 
ments of excited state BO curves [see Fig. [^. We found 
that the characteristics of LISs are extremely dependent 
on the shape of the dark state BO curve which is not di¬ 
rectly dipole coupled with the ground state. Conversely, 
the LISs are in many cases independent of the shapes of 
the bright, dipole allowed BO curves. LISs can therefore 
be used as a measure of the dark state BO curve-shape. A 
narrow single-line molecular LIS is a signature of a dark 
BO curve parallel or almost parallel to the ground state 
[see Fig.|^a)], while non-parallel curves and curves with 
shapes different from the ground state curve give horizon¬ 
tal few-line or defuse LISs. We see a continuous change 
in the molecular LISs from narrow to defuse, when more 
and more vibrational dark states are needed in the ex¬ 
pansion of the nuclear wave function of the ground state. 
The impact of nuclear motion can therefore be divided 
into characteristic categories including single-line molec¬ 
ular LISs when the expectation value of the internuclear 
separation on the dark state curve / R\Gd{Rjt)\‘^dR is 
constant in time and defuse LISs when the dark state 
curve is dissociative. 


In addition to the results shown in the present work, 
we have performed calculations on the 3D H 2 model of 
Ref. [To]. It was, however, difficult to unambiguously 
identify LISs in the spectra. There are several reasons 
for this. First, the only dark state curve (EF) in the en¬ 
ergy range of interest is very different from the ground 
state curve (X) (For information about the BO curves 
in H 2 see Ref. m)- The ground state curve is a single¬ 
minimum curve with a minimum at R = 1.37 while the 
dark state EF curve has two minima at R = 1.91 and 
R = 4.41, respectively. As a consequence the overlap be¬ 
tween vibrational ground states of the two BO curves are 
much smaller than one, and we expect the corresponding 
LISs to be very diffuse. In contrast to the ATAS spectra 
of this paper, the spectra of H 2 are very rich in structures 
(see Fig. (5) in Ref. [lO]). Two bright state BO curves (B 
and C) are located in the energy range where we expect 
to see LISs and as a consequence it is nearly impossi¬ 
ble to distinguish LISs from, e.g., which-way interference 
phenomena. To see clear LISs the dipole matrix elements 
between the ground state and a bright state and between 
this bright state and a dark state should both be large [see 
Eqs. ^ and compared to dipole matrix elements de¬ 
scribing processes responsible for other spectral features 
in the energy region of the LISs. In H 2 unambiguous 
identification of LISs is hindered by which-way interfer¬ 
ence between vibrational states of the B and C curves 
since the dipole matrix elements are almost ideal for this 
process. Tuning the coupling strengths of the system as 
outlined above, to circumvent the difficulties of H 2 , we 
observe LISs behaving very similarly to what we saw in 
Sec. |III[ We therefore expect that also in experiments 
LISs may in many cases be observed and used to divide 
nuclear dynamics on dark BO curves into categories: sta- 
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